Enhanced violation of a Leggett-Garg inequality under nonequilibrium thermal conditions 
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We investigate both analytically and numerically violations of a Leggett-Garg inequality (LGI) for a compos- 
ite quantum system in contact with two separate reservoirs at different temperatures. Remarkably we find that 
LGI violations can be enhanced when a heat current is established at low temperatures in a steady-state regime. 
Based on a Kraus operator decomposition of the non-unitary evolution for a system formed by two interacting 
spins or quantum bits, we provide analytical support for power law relations between dissipation strength and 
mean temperature in the borderlines separating parameter regions where non-equilibrium conditions affect dif- 
ferently the maximal LGI violation. Furthermore, a correspondence between spatial and temporal correlation 
inequalities is shown to persist even in such nonequilibrium thermal settings. 
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Since the seminal work by John Bell [ 1 ] about quantum 
spatial correlations, an ever increasing number of experiments 
has been devoted to test the foundations of quantum mechan- 
ics. In a parallel way, and triggered by those tests, practical 
schemes for conveying quantum information have flourished. 
New aspects of the fascinating world of quantum correlations 
arose when temporal, instead of spatial, correlations were pro- 
posed by Leggett and Garg [2], and subsequently by many 
other authors [ 3 - 5 1. In those works different inequalities be- 
tween two-time correlations have been proposed that should 
hold whenever a classical description is valid. An emerging 
Leggett-Garg inequality (LGI) violation would mark a bor- 
derline between the classical and quantum worlds. Of special 
interest has been the guide provided by LGI for testing macro- 
scopic realism in condensed matter systems l6lfT2) as well as 
in quantum optics setups 1021 [14) where focus has been put 
on LGI violations for open quantum systems in contact with 
realistic reservoirs. 

On the other hand, one of the most important challenges 
to detect superposition of quantum macroscopic states is the 
fragility of these states caused by decoherence effects. De- 
spite the ubiquitous presence of non-equilibrium situations in 
quantum mesoscopic and macroscopic physics, few formal re- 
sults about quantum correlations in such regimes are available 
1031118) . Many efforts encompassing a large variety of phys- 
ical systems have been directed to study the effects of spatial 
correlations, as measured by the concurrence [0210*?, 20 ] and 
quantum discord (QD) [21 ], under non-equilibrium gradients. 
However, LGI studies have been very scarce in systems out 
of thermodynamic equilibrium, although time correlations in 
condensed matter systems show complex and interesting be- 
haviors, like those observed in ultracold atom systems (22), 
inductively (capacitively) coupled flux (charge) superconduct- 
ing qubits [23], and even energy transfer processes in photo- 
synthetic purple bacteria 124) . 

How do thermal non-equilibrium conditions affect quantum 
two-time correlations? From a theoretical point of view im- 
portant insights have been given in Ref. 125) where nanother- 
modynamics relevance of quantum coherences has been ad- 
dressed. From the experimental side Ref.[26] shows an ele- 



gant experiment where thermal and electric currents interfere 
in a Josephson junction. Here, we propose a set up where 
the relationship between two-time quantum correlations and 
out-of-equilibrium thermal conditions can be systematically 
explored through the LGI violations. Here, we examine a LGI 
for a two-qubit system and show that violations are enhanced 
when a heat current, J, is flowing through it. 

For nonequilibrium dynamics, what is the optimal rela- 
tionship between dissipation strength (r) and temperature 
range (T) which allows the detection of quantum correlations 
through macroscopic transport measurements? The theoreti- 
cal approach we provide, in terms of Kraus operator dynamics 
evolution, gives a general framework to identify non-classical 
features through LGI violations for quantum systems in con- 
tact with two memoryless (Markovian) separate reservoirs at 
different temperatures. We find that dissipation, inter-qubit 
interactions and mean temperature compete, resulting in a 
rich T-Y phase diagram where power-law relations mark the 
boundaries between different zones where LGI violations are 
either absent, decreasing or increasing with an applied tem- 
perature gradient. 

The two-time correlation inequality, LGI, extends the well- 
known Bell inequality (designed to probe spatial quantum cor- 
relations) to the time domain. One form the LGI can take is 
(USE): 

F(t u t 2 ,t 3 ) = C(t u t 2 ) + C(t 2 ,t 3 ) - C(t u t 3 ) < 1 (1) 

where C(U,tj) is the two-time correlation of a dichotomic ob- 
servable Q (eigenvalues a=±l) between times ti and tj and 
ti < t 2 < £3. To evaluate F(ti,t 2 ,ts) in Eq.{T]), a Kraus 
operator K(t) approach has been used (28] [29). In the fol- 
lowing we limit ourselves to a steady-state regime. Due to 
this last assumption, two-time correlations depend only on the 
time difference between two instants, which implies the non- 
invasive measurement requirement in LGI settings [4 ]. The 
steady- state or long time limit of the density operator will be 
denoted as p ss = p(t — > 00). By choosing in F(ti,t 2 ,ts) all 
three times t±, t 2 and £3 in the steady-state regime and fixing 
them as t\ = 0, t 2 = r and ts = 2r, the two-time cor- 
relation function can finally be expressed as (see Ref. |29l ) 



2 



C(r) = 1 — 4p + + 4Re [g(r)], allowing us to rewrite Eq.{T]) by n(u) = \ n i(&) with index i representing each sep- 



as 



F(0, r, 2r) = 1 - 4p+ + 4Re [2#(r) - #(2r)]) < 1 (2) 

where p + = Tr |n + p ss | denotes the long time limit of the 

probability for obtaining the result a — +1, II + the corre- 
sponding projector and 
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This last term is of highest importance for the present work as 
it is the responsible for LGI violations. 

We consider a bipartite quantum system in contact with sep- 
arate thermal baths at different temperatures T\ and T2 . To be 
specific, the quantum system of interest is a two spin chain 
described by: 
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where is the energy splitting of the i-th qubit, cr^ de- 
note Pauli matrices and V describes the inter-qubit interac- 
tion strength. The baths are represented by sets of harmonic 
oscillators with Hamiltonians B4 = J2 n ^n,i^ n i&n,i where 
^ n i (a<n,i) creates (destroys) an excitation in the i-th reser- 
voir while the coupling of each qubit with its separate reser- 



voir is given by H int4 = erf ^ n g. 



where denote the system-bath coupling strengths. In the 
Born-Markov framework the system-bath couplings are es- 
sentially determined by the bath spectral density of the form 
Ji(u) = Trii(uj) where T denotes the coupling strength (taken 
identical for both reservoirs) and rii(u) = (e^ iUJ — the 
Bose-Einstein distribution of excitations in the i-th bath at in- 
verse temperature (HJ. For the sake of simplicity we limit 
ourselves to the symmetric- qubit case, i.e. e\ = 62 = e. In 
the following, energy, frequency and temperature will be mea- 
sured in units of the inter-qubit interaction strength V while 
time will be expressed in units of V~ l by letting h — Kb = 1. 
We thus set e = e/V, f = T/V, u = u/V, f = T/V, and 
f = Vr. 

The eigenstates for the qubit Hamiltonian in Eq.Q are 
given, in the a Zii basis, by: |1) = |+,+), |2) = (|+, — ) - 
I-, +))/v / 2, |3) = (|+, -> + |- +))/v / 2 and |4> = |-, -). 
The non-equilibrium thermal steady- state density operator for 
the pair of qubits turns out to be diagonal in this basis taking 
the form of a direct product as lfT5l 
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where pj s denotes the steady-state density operator of a sin- 
gle fictitious qubit of energy splitting ujj in thermal equilib- 
rium with a single reservoir at an effective temperature corre- 
sponding to the mean number of thermal excitations as given 



=1,2 

arate thermal bath. The two effective qubits have energy 
splittings given by &\ = \e — 1| and U2 = e + 1. Thus, 

Pj S = dia g{ 2 n^) ) +i ^ 2n(^)+ 1 i } imping that the stead Y- 
state spin chain density operator in Eq.([5]) does not depend on 
the system-bath coupling strength T. Of special interest for 
the discussion of results below are the coefficients C2 and C3 



in Eq.(5 ), given by C2 

/ n(d>i)+l \ / 71(0)2) 
^2n(u>i) + l J \ 



( n(Qi) \ ( n{ 
y2n(d>i) + l J \2n 



(£ 2 ) + l 



and C3 



(w 2 )+l y 

, which determine the heat current a 



v 2n(d> 2 ) + l 

cross the quantum system by the annihilation of an excitation 
from one reservoir and the subsequent creation of another ex- 
citation in the other reservoir through a flip-flop process in the 
two-spin chain. Notice that in the special case e = 1, one of 
the fictitious qubits has zero splitting, uj\ = 0, which is equiv- 
alent to having at least one bath at infinite temperature since 
71(a)) —> 00. In order to assess under which conditions the LGI 
can be violated, F(0, f , 2f ) > 1 in Eq.([2]), we pick the sin- 
gle qubit operator Q = a Zj i as the dichotomic observable to 
evaluate the LGI. The projector to the |+) eigenspace for this 
observable is n+ = |+)(+| <g> i = |1)(1| + ^ "£l=2 HH 

To gain physical insight on quantum two-time correlations 
and LGI violations for a nonequilibrium thermal spin chain, 
we first consider a simple case: once the spin chain reaches 
the nonequilibrium steady-state regime, we proceed to break 
its couplings to the heat reservoirs and let it evolve unitarily 
under the action of Hamiltonian Hq . In this case, the LGI-F 
function can be analytically evaluated, yielding to 

F(0,f,2f) = 1 + r 23 [2 cos (2f) - cos (4f) - 1] (6) 

which only involves, through r23 = C2 + C3, the density ma- 
trix elements corresponding to eigenstates |2) and |3), as ex- 
pected. Under the sole action of Hq the states |1) and |4) 
are thus decoupled from the other two states. It is then easily 
seen that F(0, f, 2f) begins at a value of 1 for f = 0. The 
LGI is violated during an initial time interval from f = to 
f = 7r/4, reaching a maximal violation when f = n/ 6 ir- 
respective of the bath temperatures. However, the maximum 
value of the LGI violation (MLGI), as defined by MLGI = 
p(Max) ^ 2f) — 1 > 0, does depend on r23 and conse- 
quently on the bath temperatures. 

Now, we restore the couplings of the spin chain with the 
baths which makes it imperative to evaluate the quantum 
open system's Kraus operators numerically (details will be 
presented elsewhere). Using this procedure, the two-time 
correlation function Cjv(t) can be safely approximated by 
a fitting correlation function of the form Cf(t) = 1 — 



^23 



1 - e~ Tvf cos(2f) 



where ^ 



£ PMC,-) + 1]. 
i=i 

The fit is good enough that \C N (f) - C F (f)\ < 1(T 5 for 
the full numerical results in the range f < 6, validating our 
analytical results described below. Thus, by using the fitting 
correlation functions, C f (t), the LGI in Eq.([2]) can be written 
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as 



T(0,f,2f) 



r 2 3 



(7) 



2e~ L VT cos (2f) - e 



-2T v f 



COS (4f ) 



< 1 



In contrast to the steady-state nonequilibrium thermal concur- 
rence (El [HI and QD ED, the LGI-T(0, f, 2f) depends ex- 
plicitly on the dissipation rate T. It is of special interest to re- 
late the LGI-T function with a measurable transport quantity 
such as the heat current flowing through the spin chain in re- 
sponse to a temperature gradient. An analytical expression for 
the heat current under a temperature gradient, AT = T\ — T 2 , 

as calculated from J (AT) = Tr jiTg£i j can be found in 

Ref. lfT5l . Although a closed equation relating F(0,f, 2f) 
with the heat flow is possible, the expression is cumbersome 
and will be skipped here, since their relationship is still better 
appreciated by looking at the figures as discussed below. 




FIG. 1: (Color online). LGI-F function at thermal equilibrium, i.e. 
AT = 0. T(0, f , 2f ) for symmetrical qubits with different splittings 
and identical bath temperatures T\ — Ti — 1.1. The LGI is violated 
whenever T(0, f, 2f) > 1. 

To illustrate the behavior of the MLGI violation under 
nonequilibrium thermal conditions we proceed to consider re- 
alistic parameters. The set of chosen parameters (in V units), 
for the results we discuss below, are in the actual experimen- 
tal range for superconductor qubits ~ 5 and T ~ 10 ~ 3 
l23l while for excitons in biological photo synthetic systems 
£i ~ 1.4 and T ~ 0.4 [24 j. First, we analyze results for the 
equilibrium case, i.e. AT = 0. As plotted in Fig. [I] violations 
of LGI are still visible up to f = 10 for an equilibrium temper- 
ature Tm = Ti = T 2 = 1.1. For weak inter-qubit couplings 
(e > 1) the violation of the LGI is small compared with the 
strong inter-qubit coupling result in the range of < f < 1. 
However, for 1 < f < 5 that behavior is reversed and the clas- 
sical non-violation regime prevails (T(0, f , 2f) < 1). More- 
over, in the weak inter-qubit coupling the LGI violation stands 
for longer times. By decreasing the qubit splitting e the LGI 
violation is becomes large for short times, as can be seen in 
Fig[T] For the special case e = 1, the LGI violation breaks 
down for any f and the classical behavior emerges, as should 
be for a quantum system in contact with an infinite tempera- 
ture bath as pointed out above. 

Now, the essential question is: can the LGI violation be 
enhanced by a temperature gradient or equivalently by a heat 



current? This is analyzed by plotting in Fig. [2] (a) the MLGI as 
a function of J and the mean temperature Tm = (Ti +T 2 )/2. 
The variation of MLGI vs. heat current (or equivalently tem- 




FIG. 2: (Color online). Composite quantum system in contact with 
two separate thermal reservoirs, (a) A 3D plot of the nrsf MLGI (for 
f = 7r/6)asa function of mean temperature and heat current for a 
two-spin chain with e = 3 and T = 0.05 is shown, (b) Slices taken 
from the 3D left panel are shown for different Tm . Notice that for 
the lowest values of Tm, MLGI increases with the application of a 
temperature gradient. Similar results are shown in the inset for higher 
mean temperatures: Orange (Tm = 2.25), magenta (Tm = 2.5), 
purple (f M = 2.75), black (f M = 3.0), gray (T M = 3.25). A 
qualitatively similar behavior can be observed for the two-spin chain 
concurrence (c) and quantum discord (d) as a function of the heat 
current. The color convention for curves in (b-main figure), (c) and 
(d) is the same. Since the values of e are the same, a symmetric curve 
is obtained for AT < 0. 

perature gradient) changes its concavity indicating that the 
enhancement of LGI violations under nonequilibrium ther- 
mal conditions is restricted to the low mean temperature sec- 
tor. This behavior can be understood as a competition of the 
coherent flip-flop spin processes (dominated by the V term) 
and a temperature dependent decay rate (TV), as can be de- 
duced from Eq. [7] Since Ty « 1, it is possible to ob- 
tain an analytical expression for the MLGI as MLGI ~ 

^23 e -f le 3^ — 1 9 that depends on the product of 
the stationary spin chain level populations r 2 3 and an effec- 
tive decay term. For low Tm, and AT << 2Tm, the popula- 
tion term in the last analytical MLGI expression is much more 
important than the effective decay indicating that the MLGI 
starts to grow as the r 2 z increases for AT = 0. By contrast, 
for higher Tm it is very easy to demonstrate that r 2 3 decreases 
faster than the contribution of the effective decay Ty indicat- 
ing that the decay term dominates and the flip-flop interaction 
does not contribute appreciably due to the fact that the popu- 
lation of each eigenstate is saturated at 1/4. This is the main 
reason why moving away from thermal equilibrium the MLGI 
decreases for high Tm • The non-equilibrium enhancement of 
spatial quantum coherences in the same spin-chain system, 
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as measured by the concurrence 1 15 ] and QD [21 ], is shown 
in Fig. |2](c)-(d), respectively. This correspondence between 
spatial and temporal correlations extends to a nonequilibrium 
thermal setting previous observations of a perfect mapping be- 
tween Bell inequality and LGI l30ll . 

Our method as applied to the single spin case (4| [6) (re- 
sults not shown here) demonstrates that MLGI is not en- 
hanced, rather it is suppressed, when a nonequilibrium tem- 
perature gradient is applied. Therefore, a two-qubit system is 
the smallest spin chain where nonequilibrium enhancement of 
LGI violations can be observed. 

A complete presentation of LGI violations under nonequi- 
librium conditions is given in Fig [3] (a) We point out that ther- 




FIG. 3: (Color online), (a) Contour plot of MLGI in the (f, T M ) 
parameter space for e = 3. The numbers around the contour lines 
give the values of the MGLI at AT = 0. Inset: (f , Tm) in log- 
log scales showing power laws for the borderlines between different 
non-equilibrium MLGI behavior regions, (b) Typical behaviors of 
concurrence and QD at equilibrium (both normalized to their maxi- 
mum value), where a blue line indicates that a temperature gradient 
increases the quantum spatial correlation while a yellow line indi- 
cates that a temperature gradient decreases the quantum spatial cor- 
relation. 

mal nonequilibrium effects on LGI violations can be grouped 
in three zones in the whole (T, Tm) parameter space. First, 
there are no LGI violations, both in equilibrium and nonequi- 
librium, for sufficiently high Y or Tm values, corresponding 
to the upper right white zone in Fig. [3] (a). This behavior 
is consistent with the deleterious effect of strong coupling or 
high temperature on quantum correlations. Second, there is a 
zone with intermediate values of Y and Tm (yellow region in 
Fig[3] (a) where MLGI > at thermal equilibrium, but the 
MLGI starts deteriorating rather than improving as the system 
is subject to a small temperature gradient. As shown in the 
inset, where a log-log scale has been used for Y and Tm, the 
numerically determined borderline between the white and yel- 
low zones has a power law expression such that Tm ~ 



This simple relation can be easily justified by using our ana- 
lytical approximation for MLGI, by searching for the (T, Tm) 
values yielding to a null value for MLGI at AT = 0. Third, 
and most importantly, there is aparameter region of low Y 
or Tm values (blue zone in Fig|3] (a) where the equilibrium 
value of MLGI is further increased as the system is taken 
out of equilibrium by setting a temperature gradient. The nu- 
merically determined borderline between the yellow and blue 
regions in Fig [3] (a) satisfies also a power law of the form 
Tm ~ T -1 / 3 . This behavior can also be explained by re- 
quiring that fgfpr l A f=o = 0, to obtain|29l 



1 

M 



(8) 



in agreement with the numerical results plotted in Fig. [3] (a). 
Clearly, when V > e no solution for Eq.([8]) exists indicating 
that a requirement for the thermal nonequilibrium enhance- 
ment of LGI violations is that spins are in the weak coupling 
regime. By contrast with two-time quantum correlations, two- 
point spatial correlations in the steady- state do not depend on 
the quantum system coupling strength with the thermal baths. 
However, as shown in Fig|3](b), the two- spin concurrence and 
QD also present a nonequilibrium thermal enhancement at 
low enough temperatures signaling one more time the corre- 
spondence between temporal and spatial quantum correlation 
behaviors. As it has already been documented the QD persists 
up to larger temperatures as compared with the concurrence. 

In summary, we have provided a microscopic derivation 
of LGI, for a spin chain coupled to two reservoirs at dif- 
ferent temperatures. Nonequilibrium thermal quantum cor- 
relations have been analyzed exactly by mapping the origi- 
nal two-interacting spins coupled chain to an equivalent sys- 
tem of two-noninteracting spins, each one of them in contact 
with a renormalized heat reservoir at an effective temperature. 
Based on a Kraus-operator approach, we demonstrate that in 
a certain range of temperature gradients the steady- state LGI 
violation can be enhanced. The frontiers between different 
behaviors of the MLGI response to thermal nonequilibrium 
conditions have been found to be characterized by power laws 
relating mean temperature of the heat reservoirs and the spin 
chain coupling strength with the baths. This leads to the inter- 
esting feature that nonequilibrium thermal conditions provide 
the opportunity to enhance not only spatial but also temporal 
quantum correlations. 
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